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1.  Introduction 

In  the  final  report,  given  no  consideration  to  the  modern  achievements  in  the 
problem  under  study  and  to  the  analysis  of  a  number  of  physical  processes  to  justify 
the  problem  statement  that  are  given  in  detail  in  reports  of  Items  2-4,  we  present  only 
those  sections  in  which  the  formulae  are  directly  used  in  two-dimensional  code.  The 
algorithm  structure,  program  code  and  its  description,  and  also  some  results  of  cal¬ 
culations  are  presented  in  appendix. 


2.  Calculation  Methods  of  Electron  and  Ion  Temperatures  and  Charge 

Composition 

Since  in  kinetic  equations  for  calculation  of  charge  composition  and  in  the  equation  for 
electron  temperature  the  terms  connected  with  non-elastic  processes  have  the  same  structure, 
clearly  the  differential  algorithm  should  be  built  from  the  integration  of  blocks  that  calculate 
their.  In  three-speed  approximation  the  system  of  ionization-temperature  equations  takes  the 
form: 


+div(nlz*Vl)  =  Slz 

(2.1) 

rh 

+ div(n2z  *V2)  =  S2z 
a 

(2.2) 

Z=0,l,2,...Zm 

/fo 

-^  +  div(ne*Ve)  =  Se 

(2.3) 

0  rtrp 

-*«1  *rl  *div^  =  01 

(24) 

*n2*  2  +  n2*T2*divV2=Q2 

2  at 

(2.5) 

— ~  +  *T  * divV  -  div(Z,VT,)  +  S„  +  Q, 
at 

(2.6) 

Calculation  has  shown  that  ion  contribution  to  variation  of  electron  density  is  infini¬ 
tesimally  small  in  comparison  with  electrons  themselves.  Therefore,  the  velocities  SiZ  for  in- 


coming  flow  and  S2Z  for  reflected  flow  are  defined  by  collisions  of  ions  with  electrons  and  by 
the  overcharging  processes  between  ions  of  counter  flows. 


Slz  =  *ne*  jz_Xe  - nlz  * n)  *  jez) - (nXz  * ne  *  jze  - nXz+x  * n)  *  jez+ , ) - 
-  «iz  *  ne  *  fez  +  nu+]  *  ne  *  jvez+x  -  nu  *  n2z+x  *  j\2  -  nXz  *  n2z_x  *  j£  + 

+  n2z  *nlz_x  *j2X  +n2z  *nXz+x  *  jz2X 


(2.7). 


The  similar  expression  is  for  S2z.  A  detail  analysis  of  reaction  rates,  expression  for 
variation  rate  of  electron  temperature  in  non-elastic  processes  See  and  expressions  for  energy 
transfer  in  elastic  collisions  Qi,  Q2  and  Qe  are  presented  in  the  previous  report.  It  is  signifi¬ 
cant  that  the  variation  of  Te  due  to  yield  of  photorecombination  radiation  is  already  taken  into 
account  with  sufficient  accuracy.  The  terms  with  braking  and  magnetic  braking  radiation  are 
added  in  See  by  separate  terms,  since  they  did  not  connect  with  non-elastic  kinetic  processes. 

Electron  thermal  conductivity  is  defined  by  the  coefficient: 

,  (■ k% ,)M 


m 


\£*Z*e**  In  A 


;  i/cm*c. 


(2.8). 


From  the  simple  estimations  it  follows,  that  the  characteristic  distance  of  heat  propagation  by 
electron  thermal  conductivity  in  a  time  of  the  interaction  of  incoming  and  reflected  flows 

2  A,, 

(x~4*icr8  s)  is  defined  by  a  coefficient  of  temperature  conductivity  =  —  *  —  and  has  the 

3  n. 


value  L  =  ^J%e  *  t  comparable  to  CT  radius.  However,  taking  into  account  that  the  scheme 
itself  has  a  differential  dissipation  and  in  our  attempt  to  understand  in  more  detail  the  features 
of  interaction  in  a  process  of  the  interpenetration  of  incoming  and  reflected  flows,  we  are  re¬ 
stricted  only  by  the  effect  of  differential  dissipation. 

The  basic  requirements  that  we  posed  at  development  of  the  kinetics-temperature  block 
are  autonomy,  universality  and  economy  of  calculation  in  time.  The  autonomy  should  allow  to 
use  the  block  for  solving  of  extended  problem  categories,  essentially  independent  of  gas  dy¬ 
namics  flow  pattern  and  of  the  chosen  grid  -  type  of  Lagrangian  or  Euler.  In  this  case  the  par¬ 
ticle  velocities  and  densities  should  be  by  input  parameters  for  calculation  (  aiz,  Te,i )  at  a  fu¬ 
ture  time  layer.  The  block  universality  should  provide  a  possibilty  of  calculation  (diz,  Te,i  ) 
essentially  in  any  wide,  but  reasonable  range  of  the  variation  of  initial  data:  n<!022  cm'3,  T= 


0.025  100  KeV.  Since  the  implicit  differential  schemes  are  on  the  base  of  calculation  algo¬ 

rithm  (  aiz,  Te,i ),  then  the  special  attention  is  focused  on  the  convergence  of  iteration  process  - 
its  rapidity  and  economy  in  time. 

Since  the  equation  for  aiZ=niZ/ni  and  a2z=n2z/n2  are  linear,  in  the  general  case  of  the 
no  quasineutrality  the  equation  for  a=ne/(nx+n2)  is  calculated  separately  and  the  iteration 

;2z  a  is  created  that  is  rapidly  converged  under  ni,  n2  < 

5.*10x9  cm-3.  In  our  problem,  as  shown  in  numerical  experiment,  some  deviation  from  the 
quasineutrality  is  observed  only  near  sharp  fronts.  In  general  the  quasineutral  approximation  is 
applicable  and  the  iteration  are  carried  out  with  the  relation: 

a  =  YJZ*au+YJZ*a2l  (2.9). 

z  Z 

However,  the  iterations  are  not  satisfied  to  universality  criteria,  since  at  n>  s.no^-r-io20  cnr3 
their  convergence  slows  down  sharply  at  T<  1  eV.  Therefore  the  special  iteration  processes 
was  developed,  which  in  the  field  of  single  ionization  at  high  particle  density  up  to  1022  cnr3 
produced  the  rapid  convergence.  So,  the  kinetic  calculation  block  (aiz,  a2z,  a)  in  the  last  ver¬ 
sion  is  satisfied  to  all  requirements  outlined  above. 

The  similar  requirements  are  on  the  base  of  the  production  of  the  block  of  temperatures 
Te,  T„  T2  calculation.  We  proposed  an  idea  of  decomposition  of  temperature  equation  on 
physical  processes.  Through  introduction  of  so  called  temperatures  yi,  y2,  ye: 


yi-Ti/Tig,  y2-T2/T2g,  ye-Te/Teg  (2.I0), 

in  temperature  equations  we  can  decompose  the  terms  responsible  for  adiabatic  expansion 
(«,  *Ti*divVj )  from  ones  corresponding  to  non-elastic  and  elastic  collisions  (Qi,Q2,Qe,See). 
Substitution  of  (2.10)  in  (2.4-2.6)  gives  for  adiabatic  temperatures  Tig,T2g  the  equations: 

3.  dTiv  - 

—  *  ~^L  +  Tig  *  divVi  =  0. ,  r«e  1=1,2,  (2. 1 1), 


that  in  combination  with  continuity  equations: 


^3/2 

have  the  solutions  — —  =  const.  In  doing  so,  for  kinetic  temperatures  yi  the  equations  are 

ni 


taken  without  divergent  terms: 


dy,  g 

dt  3  n,  *  Tjg  ’ 


TJSfi  i=l,2 


(2.13)- 


The  continuity  equation  for  electrons  (2.3)  has  a  non-zero  right  hand,  therefore  substituting 
(2.10)  in  (2.6)  we  obtain  the  equation: 


3  dv 

—  *n  *T  *-^-  +  n  *T  * 
2  e  eg  dt  e  e 


3  1  dTee  l  dn, 

—  *  . 


v2  Teg  dt  ne  dt 


+  Te*Se  —  Qe  +  Se 


(2. 14). 


Equating  the  expression  in  parenthesis  to  zero,  we  get  the  decomposed  equation  system  for 
calculation  of  electron  temperature: 


§*",*r„.%=e,-7;.s,+s„ 

•  T,n  (215). 

=  const 

.  He 

The  described  idea  of  decomposition  is  convenient  to  the  Lagrangian  grids,  since  it 
allows  one  to  transform  the  equation  system  in  partial  derivatives  to  a  system  of  regular  dif¬ 
ferential  equations.  However,  with  the  Euler  grid,  using  the  “method  of  big  particles”  as  it 
adequately  describes  physics  of  gas  dynamics  processes  ,  the  Lagrangian  type  of  calculation 
corresponds  to  the  described  decomposition  method.  Without  introduction  of  decomposition 


method  the  terms  («,  *T;  *divVj)  introduce  additional  approximation  errors,  since  the  used 
form  of  temperature  equations  (2.4)-(2.6)  is  not  a  divergant-conservative. 

It  should  be  emphasized  that  main  calculation  difficulties  were  connected  with  the 
terms  Qi,Q2,Qe,See  in  temperature  equations.  This  needed  the  development  of  implicit  differ¬ 
ential  scheme  and  organization  of  special  iteration  process,  this  work  is  completely  performed 
in  the  frame  of  given  contract.  In  calculations  the  main  result  was  achieved:  rapid  conver¬ 
gence  of  iteration  process  at  reasonable  accuracy  (~  1%). 


3.  Momentum  and  Energy  Transfer  Rate 

For  a  wide  range  of  plasma  current  conditions  a  local  Maxwellian  particle 
velocities  distribution  is  applicable  in  terms  of  their  directed  motion  velocities.  As  a 

result  of  the  detail  analysis  we  obtained  the  rather  complete  expressions  for  Ri2  and 


Q  some  of  which  were  already  used  to  solve  the  problems  of  clot  collisions.  With  these 
expressions  for  further  studies,  we  show  them  here  in  the  rather  complete  form: 

^  =  Ej,  +  R,o  +  ^  (3-1) 

First  term  corresponds  to  collisions  with  ions,  second  -  with  neutrals,  and  third 
-  with  electrons: 


-  Airy  4  *  72  *  72 
p*lF 

R,c=~ 

V71 

(3-2), 


K  =  — 4  * \ 1* *1/ *ne  *(K~K)*Fe K ) 
«  *if? 


(3-3)- 


The  expression  for  has  a  similar  form.  For  electrons  we  get: 


R.  =~ 


4ne4L  Z?*nu 


mB 


*nt  *[- 


W? 


HK-V,)*Fe(xu)  +  ^^*(K-V2)*Fe(x2e)] 


if, 


It 


(34). 


Here: 
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IF 


■yj2kTx  !mx  +2^7^  /#?, 


:;«  =  ^2^7^  Jmx  +~2kT2  jm2\ 


X,.  =K  /p kTJm,iX2,  =r2,  IpkX  bn. 


(3-5). 


The  ct12  is  the  ion-neutral  collisional  cross  section.  For  energy  transfer  rate  we  get: 
Qi=Qii+Qio+Qie ,  where: 


4to4Z2Z2L  |i(T2  -Tx)  x 

a,  = - —^nXin2i  *  - —  *  <D(x)  + 

p*lF 

^27K  +  WiT2K  W  2x  2 mi 

2 1  '  i-"^*[°(x)-^exP(-y)]} 

m2Tx+mxL2  W  yjn 


j2io  - 


JSs'l,..,.  ,,KT2-t,)  r,K  +  »,T2K.ir 


»■  *  »20  *  <  ^  V'l'Wt  z  1 '  'T 

mxl2+m2Tx  m2lx+mxL2  a 


*-r*FoW) 


<2„  =  **„*».*  r-(T>  -4  *  ®(*„ )  +  » [®(xi, )  -  *  exp(-4 ))} 


(3-6). 


Here:  W=V2-Vt;  W,,  =V, -V,;W  =  \w\;  W?,,  =\wu  . 


The  expression  for  Q2has  a  similar  form.  For  Qe  we  find  (with  me«mi ): 


k=\  Wke  mk  Te 


V.  *W.e 


(37). 


Here  the  summation  is  over  ion  kind  k=i,2.  The  dimensionless  function  arguments  in 
right  parts  and  their  asymptotic  forms  look  like: 


2  * 

0(x)=  — ■==*  [exp(-/2)^/=< 
Vit  o 


—j=  *  (x- x3  /3  + . );  x-»  0 

■yjn 

.  exp(-x2) 

1 - h= - -  + . ;  x-»  oo 


[  Vit  *  x 

r  4 

-=  *  (x3  /3  -  x5  /5  + . );  x->  0 

yjn 

[l  -  x*  exp(-x2)  /4tz  + . ;  x— »  oo 

r—  11  1  o  f8  /3  +  8x2  /15;  x— >  0 

F0(x)  =  y/n*x*(l O(x)  +  (1  +  — j)  *  exp(-x  )  =  \ 

Ax  x  2x  [yln*x,x-J>  oo 


F,-  ( x )  =  Fe  (x)  =  ®(x)  — =  *  exp(-x2)  =  <j  yfn 


'F(x)  =  yfn  *  ( x3  +  3x+  —)  *  ®(x)  +  (x2  +  5  /2)  *  exp(-x2)  = 

4 


4  +  9  *  x2;  x— >•  0 

Vit  *  x3;x— »oo 


(3-8). 


The  analysis  of  cross-sections  and  rate  coefficients  is  presented  in  the  Item  2  of  the 
report. 

The  system  of  the  dimensionless  kinetic  equations  looks  like: 


dcq 

dt 


—  A00a0  +  A01oij; 


da,  —  _  —  _  —  _ 

—  A10a0  +  AjjOj  +  A12a2; 
dt 


^am  -  a  &  +  a  a 


(3-9). 


It  is  a  similar  system  for  az.  The  initial  structure  was  determined  on  the 

quasibalance  model.  The  values  of  coefficients  Ay ,  Ay  (one  line  corresponds  to  “l”, 

two  line  correspond  to  “2”)  are  determined  by  overcharging  rate  and  electron  kinetics. 
Then  calculations  from  the  kinetic  block  determined  the  ionization  degree  and 
average  charges: 


n,  -A  n,  ^  —  n.Z.  +  n,Z, 

ae  =  -U£zaz  +  -^*£Zaz  =  -LJ - — ; 

n  1  n  V 

m  m  __ 

n  =  ni  +n2;Zj  =  £Z*az;Z2  =  ^Zaz 


n 


(3-io). 


m  m  

During  the  computation  the  conditions  ^ot;  =  1 ;  otj  =  1  are  verified;  the 

i  i 

error  did  not  exceed  o.i%.  The  kinetic  block  for  calculations  of  ion  level  populations  is 
a  similar  kind,  determining  line  spectra  of  plasma  radiation. 


4.  Calculation  of  Linear  Radiation  Method 

the  density  mzk  of  the  ion  component  of  type  i  (i=i,2)  and  charge  of  Z  that  is  excited  to 
the  electron  state  level  k  is  determined  by  the  equation: 


%-  +  div("a 


♦ U,)  =  Sa 


(4-i)- 

Here  Sizk  is  the  rate  of  the  density  mzk  change  in  inelastic  processes: 

sizk  =  ne  * ( *  j2  km  - nizk  *  XA>,)  +  (n2e  * «,>+i  *  Je,,+i,k  ~ne*  nizk  *  jzM)  + 

m*k  m*k 

•+  ne  +  *nizn,  ~nizk  *  YaA«,kn, 

m>k  m<k 

(4.2). 

The  first  term  in  brackets  corresponds  to  transitions  between  excited  states  as 
a  result  of  electron  collisions;  the  second  term  corresponds  to  a  capture  of  free  elec¬ 
tron  by  ion  (Z+i)  to  the  level  k  as  a  result  of  triple  collisions  and  the  reverse  process  of 
ionization  by  electrons;  the  third  term  corresponds  to  the  process  of  photo-  and  die- 
electron  recombination  to  the  level  k  of  ion  (Z+i);  last  two  terms  describe  the  change 
of  n^k  as  a  result  of  spontaneous  transitions  between  electron  state  levels.  The  esti¬ 
mations  show  that  for  transitions  between  excited  states  the  plasma  is  optically  trans¬ 
parent,  while  for  transitions  to  the  basic  sate  it  is  necessary  to  take  into  account  the 
radiation  reabsorption.  The  population  of  lower  excited  state  has  an  especially  strong 
effect  on  populations  of  excited  states.  The  transition  from  this  state  (2-»i)  corre¬ 
sponds  to  so-called  resonance  radiation.  Usually,  to  take  into  account  a  self¬ 
absorption  of  the  resonance  radiation,  the  approximation  of  Biberman-Holstein  [2,3] 
is  used  by  the  inputting  of  photon  escape  probability  ©z  from  a  plasma  volume. 


We  obtained  the  sufficiently  general  expressions  for  ©  by  the  example  of  ex¬ 
panding  plasma  volume,  when  the  plasma  layers  moved  relative  to  each  other,  and 
sufficiently  convenient  for  calculations  the  approximation  function  is  given: 
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A  3 
A  +  —  *T 
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(A  3  1 

A  +  —  *T 

7  J 


3/7  *t 
( 


V3+, 


V 


(4.3). 


Here  a  =  2  *  0)0  *  is  the  ratio  of  the  Doppler  shift  to  the  line  half-width;  a>0  is 

the  frequency  at  the  line  center;  U0  is  the  relative  speed  of  plasma  flows;  c  is  the  light 


speed;  A co  is  the  effective  line  width;  R  is  the  effective  radius  of  plasma  volume  V 

I3  *  J/ 

(R~  3 - );  r  is  the  effective  distance  from  plasma  center  up  to  point  of  radiation. 

V  4  *7T 

In  particular,  from  (3)  it  may  be  concluded  that  plasma  motion  begins  to  effect  on  © 
at  A  >  yfr  .  At  r  =  0  we  have  ©  =  1  and  at  r=R,  that  is  on  the  plasma  boundary  (at 
A  «  t  » 1 ),  we  can  obtain  ©  =  .  So,  slightly  more  than  the  half  of  radiation 


freely  yields  the  plasma  clot  region.  At  A  ->  0  (that  is  at  U0  0)  the  expression  (3) 


transforms  to  known  expression  for  plasma  at  rest,  but  takes  into  account  the  dis¬ 
tance  to  plasma  volume  boundary  more  precisely.  The  expression  for  t  correspond¬ 
ing  to  transition  k-»i  of  ion  Z  takes  the  form: 

r  =  — —  *  ^zk'  *n  * R 
4  Ao)z  2 

(4-4). 

where  Xz  kx  =  1 .24  *  1 0"4  /  Ezk  (eV) ,  cm  is  the  wavelength  of  transition  k-»i  for  ion  Z;  Ezk 


is  the  energy  of  level  k  relative  the  basic  state. 

The  expressions  for  probabilities  of  spontaneous  and  collisional  transitions 
were  presented  in  the  Item  4  of  the  report. 

All  plasma  optical  characteristics  and  first  the  intensity  of  its  linear  radiation 
may  be  determined  by  solving  a  system  of  linear  equations  relative  to  aizk  for  each  Z: 


(4-5). 


««*  * 
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£  JzM  +  JzM  +  — —  *  £  AzM  *  0 
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10  /  1  > 

-  £  atm  *  Am  +  — —  *  Am  (m  >  *) 
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=  «,,z+i  *(oc*n*  jez+xk  +  jve2+lk ) 


The  power  radiation  in  the  given  line  is  defined  by  expression: 

Sum  =  Am  *  Ezm  *  ntk  >  eV/cm3*s. 

In  the  region  of  the  existence  of  incoming  and  reflected  flows: 

=  Am  *  EzM  *  (nizk  +  «2  zk  )  =  Am  *  EzM  *  (aizk  *ni+  a2zk  *«2>  = 

=  Am  *  EzM  *  azk  *(«i+«2)  =  Am  *  Ezm  *azk*n,  eV/cm3*s. 

In  these  transformations  it  is  supposed  that  alzk  =  a2zk  =  azk ,  since  they  are 
determined  primarily  by  electron  temperature  and  density.  The  slight  difference  in 
aUk  and  a2zk  obtained  in  calculations,  plays  no  important  role  in  the  analysis  of 
plasma  linear  radiation. 


5.  Statement  of  Dynamic  and  Electromagnetic  Problem 
5.1.  Motion  and  temperature  equations 

It  is  supposed  that  ions  of  all  charges  of  incoming  component  have  velocity  Vx , 
ions  of  all  charges  of  reflected  component  -  velocity  V2 ,  electrons  -  velocity  Ve .  Then 
for  ions  of  charge  Z  of  incoming  component  the  motion  equation  takes  the  form: 


mx  *nz  * 


f  dVx 
dt 


+  (VXV)VX 


-  1  -  - 
=  -VP,  +nz  *Zx*e*  E  +  -*VxxB 
1  1  V  c 


+  Rn 


(5-1)- 

where  Rz<  =  £PZ[  Zj  +RZuC  is  the  friction  force  of  ions  Zi  with  all  ions  of  the  second 


component  and  electrons.  This  value  is  proportional  to  RZuZl  ~  Zx  *  niZ  *  Z\  *  n2Z . 

Since  all  ions  with  Zi  have  velocity  Vx ,  we  can  obtain  the  sum  in  friction  force  by 
summing  equations  (5.1)  over  all  Zi. 


Zz>2*”z,  HEZ2*"z2  =  rc,  * «2  * (XZ12  * «z )* (Zz2  * «Z )>  where 


Vz. 


V  £i 


^  "z 


az  = 
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_  fl 

az  =  — ,  and  in  the  electrodynamic  term  you  can  find  the  linear  sum 


£z,*nZi  =  «,  *  £z,  *az  .  The  calculations  showed  that  the  difference  under  root- 

Z\ 

men-square  and  linear  averaging  is  negligibly  small: 
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Because  of  this  the  z  is  used  in  all  calculations. 

In  such  a  manner  the  dynamic  equations  take  the  form: 
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\ 
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(5-2) 


(5-3). 


Here  i=i,2.  We  have  mi=m2=m;  Pi=ni*Ti;  Pe=ne*Te  for  incoming  and  reflected  flows. 


The  full  expressions  for  Rt,Re  and  also  the  equations  for  Ti,  T2,  Te  and  their  terms  are 
given  discussed  above. 


5.2.  Electromagnetic  field  equations 

The  complete  equation  system  describing  the  dynamics  of  three-fluid  plasma 
may  be  obtained  by  addition  of  Maxwell  equation  system  to  dynamics,  temperature 
and  continuity  equations: 

dB  * 
dt 

(54) 

dE  .  - 

—  =  — 4*n *  1  +c*  rotB 
dt 


(5.5) 


divE  =  4  *  n  *  pe 


(5-6) 


divB  =  0. 


(57), 

where  the  densities  of  electric  current  j  and  charge  pe  are  defined  by  plasma  motion 


j  =  Vl*e*'£jnZi  *Zx+V1*e*YJnz1  *Z2-e*ne*Ve  = 
=  e*nl*Zl*Vl+e*n2*Z2*V2-e*ne*Ve 


(5-8) 


pe  =  e  *  nx  *  Z,  +  e  *  n2  *  Z2  -  e  *  ne 


(5-9)- 

When  large-scale  low-frequency  plasma  motions  are  considered  rather  than  internal 
plasma  vibrations,  the  electron  generation  play  no  essential  role  in  formation  of 

(E,  B)  fields,  since  the  plasma,  by  any  way,  tends  to  be  a  quasineutral.  Therefore  in 


typical  plasmodynamics  problems  the  field  E  is  determined  from  motion  equation 
for  electrons,  and  then  the  charge  density  pe  (if  it  is  necessary  to  estimate  pe  for 


sharp  fronts)  is  defined  from  the  Poison  equation  (5.6).  By  doing  so,  the  quasineu¬ 
trality  condition  may  be  used  in  dynamics  and  kinetics  equations  with  high  accuracy 

while  divE  *  0. .  This  is  so-called  plasma  approximation  that  we  use  in  the  given 
problem. 


In  the  region  of  intensive  interaction  where 


V  -  V 

r  i  r  e 


<  1 . ,  we  can  mark  the  main 


dependence  of  velocities  in  the  expression  for  Re 


R  = —m  * n  * 

e  e  e 


rKzvL+K-v2} 


V 


'el  J 


(5-io), 

where  reX,zel  were  taken  from  the  general  expression  for  Re  obtained  and  given 
above.  Such  presentation  is  required  only  to  obtain  the  complete  equation  system  for 
field;  the  strict  expression  for  Re  was  used  in  the  dynamics  equation.  In  such  a  man¬ 
ner,  we  can  obtain  from  (5.3) 


1  -  -  VP 

E  =  --*[VexB ] - ^ 

c  e*n„ 


^Kz5  i  Yizli 

\  Tel  T  e2  J 


Substituting  Ve  from  (5.8)  in  (5.11)  we  obtain: 


£  =  --*[FxP]  + - - — 

c  e*c*n , 


-  -  VP  / 

-  *  [7  x  B] - e-  +  ^-  + 

e*n„  cr 


+  e* 


f  Z 2  *  n2  Z,  *  nx 


cr, 


■(k-v,) 


(5.12). 

where  the  last  term  may  be  simplified  as  a  result  of  using  of  the  quasineutrality  con¬ 
dition  which,  as  seen  from  the  structure  of  this  term  itself,  play  no  fundamental  role. 

Here  V  =  (z,  *ni  *Vl  +  Z2  *n2  *V2)/ne  is  the  mean-charge  ion  velocity; 
crl=e2  *ne*rel/me;  cr2  =e2  *ne*te2/me  are  the  partial  conductivity,  and 
<7  =  cr,  *<r2  /(cr,  +<r2)  is  the  total  plasma  conductivity.  Expression  (5.12)  is  essentially 
the  Ohm’s  law  taking  into  account  the  Haul  effect,  in  which  the  last  term  is  charac¬ 
teristic  for  two  ion  flows. 

Neglecting  of  displacement  current  in  (5.5),  since  IfJJeJJfJ  «  c ,  we  obtain 


7  =  —  *  rotB 
4  n 

(5.13). 

Substituting  in  (5.12),  taking  into  account  (5.4),  we  obtain  the  system  of  two  equa¬ 
tions  for  determination  of  fields  (E,  B) : 

1  -  _  1  _  _  c 

E  =  —  *  [F  x  P]  + - *  [rotB  xP]  + - *  rotB  - 

c  4nen„  4  ncr 
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■  +  e* 
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dB  - 

—  =  -c  *  rotE 
dt 


(5-14); 


(5-15). 


When  substituting  the  expression  for  E  from  (5.14)  in  (5.15),  we  can  obtain  equation 
for  B  in  diffusion  form,  with  two  last  terms  in  (5.14)  that  will  be  responsible  for  mag¬ 
netic  field  generation. 

A  full  numerical  problem  solution  on  determination  of  field  B  from  equation 
(5-15)  on  dynamic  plasma  characteristic  represents  problem  for  itself  and  demands 
the  development  of  special  two-dimensional  algorithm.  In  the  present  time  the  basics 
for  such  algorithm  were  developed  with  decomposition  method  on  physical  processes. 

In  frame  of  the  present  work  it  is  planned,  in  accordance  with  contract  con¬ 
tents,  to  carry  out  the  approximate  field  calculation  with  a  number  of  simplifications 
that  are  typical  for  given  problem. 

Since  the  thermalisation  of  directed  electron  energy  is  much  faster  than  that  of 
ions  {ree  =  z„  *  Jm/M) ,  electron  gas  behave  as  continuous  medium  behind  the  re¬ 
flected  wave  in  the  absence  of  friction  forces  over  ions  and  have  Ve  =  0.  On  the  other 


hand,  if  the  friction  force  between  electrons  and  ions  is  dominant  in  equation  (5.13), 


and  the  force  between  ions  is  not  too  high,  then  «, 


*2’ 


Vi 


V, 


El 
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re2 ,  and 


_  .  m*n* (Ve  -V.)  me*ne*(Ve+V2)  ,.  .  T.  Tr  „  . 

from  equation  — - — - — — —  — 2 — - — — —  »  0  we  obtain  Ve  »  V,  -  V2  »  0 . 


'  el 


Thus  the  friction  over  ions  «i»  and  «2»  acts  in  opposite  directions  and  again  Ve  is 


close  to  zero.  It  immediately  follows  that  electric  field  due  to  low  velocities  values  Ve , 

plays  no  important  role  in  the  region  of  interpenetrating  of  incoming  and  reflected 
flows.  Magnetic  field  can  not  significantly  effect  on  plasma  parameters  due  to  short 
interaction  time.  By  doing  so,  in  the  first  approximation  the  field  problem  may  be 
considered  as  not  self-consistent,  that  is,  to  take  no  account  the  reverse  effect  of  field 
on  plasma.  Then  from  kinetics  and  dynamics  equations  (5.2),  (5.3)  (without  consid¬ 
ering  the  fields)  we  obtain  current  value  (5.8)  and  from  equation  (5.13),  that  presents 
simply  the  Biot-Savart  law,  we  obtain  field  value  along  zone  of  flows  interaction.  The 
calculations  of  external  magnetic  field,  carried  by  such  a  manner,  are  presented  in 
figures  and  in  appendix  and  represent  the  upper  estimation  for  its  amplitude. 


5.3.  Boundary  conditions 


The  initial  moment  t=o  corresponded  to  contact  moment  of  incoming  flows 
with  a  target.  Flow  characteristics  at  this  moment  of  time  take  the  form  in  dimen¬ 
sionless  units: 

a)n,=l;  0<x<a;0<r<R, 
n2=  0;  0<x<a;0<r<R, 
n2=l;  x  =  0;  0  <r  <R, 

b)  Ux  =  -1;  0  <  x  <  a ;  0  <  r  <  R , 

F,  =  0;  0<x<a;0<r<R, 

U2=  1;  x  =  0;0<r<i?, 

V2=0;  x  =  0;  0  <r<R, 

c)  Ue=- 1;  0  <  x  <  a ;  0  <  r  <  R , 

Ve  =0;  0 < jc < a;  0 <r<R, 
d)  7;  =20  eV,  0<x<a;  0<r  <R, 

T2  =20  eV,  x  =  0;  0<r  <R, 

Te  =  20  eV,  0  <x<a  \  0  <r  <  R. 

In  all  free  boundaries  P{=P2=Pe  =  0 . 

In  the  target  surface  (x=o): 

U2  (x  =  0)  =  -Ul  (x  =  0) , 

Ue(x  =  0)  =  Ve(x  =  0)  =  0, 

V2(x  =  0)  =  Vi(x  =  0). 

By  this  means,  in  frame  of  contract  all  physical  processes  connected  with  dy¬ 
namics  and  interaction  with  a  target  the  plasma  clot  oriented  to  facility  “Marauder” 
were  fully  analyzed.  The  complete  equation  system  was  formulated  and  the  two- 
dimensional  numerical  algorithm  was  developed  which  allows  one  to  determine  all 
basic  characteristics  of  the  interaction  zone  of  plasma  clot  with  target:  charge  compo¬ 
sition,  ionization  degree,  temperatures  Tx,T2,Te,  velocities  Vi,V2,Ve  and  concentra¬ 
tions  nvn2,ne  of  plasma  flows,  to  calculate  continuum  and  linear  plasma  radiation,  to 
estimate  electromagnetic  fields. 

The  program  text,  user’s  guide  and  examples  of  calculations  are  given  in  the 
report  appendix. 


6.  Conclusions 


l.  The  given  work  gives  sufficiently  full  presentation  that  is  going  and  may  be 
achieved  for  high  temperatures  under  plasma  conditions  using  the  scheme: 


Taking  into  account  the  generator  possibilities  it  is  unlikely  that  the  temperatures  are 
not  higher  than  ~io  KeV.  In  addition,  a  number  of  physical  restrictions  exists. 

2.  In  our  opinion,  the  most  interesting  is  the  interaction  region  of  incoming 
and  reflected  flows  as  a  source  of  x-ray  radiation  that  is,  due  to  its  high  penetration, 
can  effect  on  the  space  behind  target,  change  the  conductivity  and  destroy  optical- 
electronic  elements.  The  given  work  besides  all  questions  may  solve  the  problem  on 
the  source  of  such  radiation. 

3.  The  impulse  of  electromagnetic  field  in  the  meter  and  decameter  range  pres¬ 
ents  also  the  essential  applied  interest.  The  pulse  can  produce  the  interference  in 
electrical-technical  systems  operation.  However,  the  calculation  of  such  electromag¬ 
netic  pulse  is  the  complicated  two-dimensional  problem.  In  the  given  report  on  the 
base  of  approximate  numerical  solution  it  is  shown  that  amplitude  of  magnetic  com¬ 
ponent  of  such  impulse  maybe  up  to  ios-io6  Gauss. 

In  such  a  manner,  the  electromagnetic  radiation  in  different  ranges  may  be  of 
essential  interest  in  this  problem. 


Fig.l.  Spatial  and  energetic  characteristics  of  plasma  clots  colliding  with  a  target: 
To=20-30  eV,  U0=io7-io8  cm/s,  n0=iol6-io19  cm-3 
a=3  cm,  R=3  cm,  Ri=o.25  *R,  R2=o.75*R 


Fig.2.  General  structure  of  calculation  grid.  Xmax 

Values  calculated  at  time  layer  n: 

Ux(i,j),  V»(ij),  U2(i,j),  V2(i,j),  Ue(i,j),  Ve(ij) 

Ti(i,j),  T2(i,j),  Te(i,j),ni(i,j),  n2(i,j),  niZk(i,j),  n2zk(i,j),  Ezkm(i,j) 


x 


GENERAL  BLOCK  SCHEME  OF  THE  PROGRAM 


Fig-3 


Preliminary  calculation 
of  dynamic  characteristic 

(I  -  stage  of  «particle-in-cell  method»),  - 

taking  into  account  the  flows  interaction, 
calculation  of  charge  composition 
(SUBROUTINE  STAR) 

- 1  internal 

iterations 


Determination  of  substance  flows 
through  boundaries  of  calculation  cells 
(II  -  stage  «PIC») 

Final  calculation  oflhe  characteristics" 
of  plasma  flows 
(III  -  stage  «PIC») 

- t11+1=t“  +At - 


Formation  of  files 
for  graph  drawing 


Calculation  of  populations 
levels  and  radiation 


Ult  u2,  ue, 

Ti,  T2,  Te, 

Hi,  II2,  CCZ  ,  Ilizk,  Il2zk,  Ezkm 

(Ui.dat,  Ti.dat,  ANi.dat  etc.) 


Figure  captions: 


Fig.  4  Spatial  distribution  of  ion  temperatures  in  the  incoming  and  reflected  plasma 
flows  at  different  time  moments  (one-dimensional  version  n0  =1017  cnr3). 

Fig.  5  Spatial  distribution  of  ion  velocities  in  the  incoming  and  reflected  plasma  flows 
at  different  time  moments  (one-dimensional  version  n0  =io17  cm-3). 

Fig.  6  Spatial  distribution  of  ion  densities  in  the  incoming  and  reflected  plasma  flows 
at  different  time  moments  (a)  and  ionization  degree  (b)  (one-dimensional  version  n0  =1017 
cm-3). 

Fig.  7  Spatial  distribution  of  ion  temperatures  in  the  incoming  and  reflected  plasma 
flows  at  different  time  moments  (one-dimensional  version  n0  =1018  cm-3). 

Fig.  8  Spatial  distribution  of  ion  velocities  in  the  incoming  and  reflected  plasma  flows 
at  different  time  moments  (one-dimensional  version  n0  =lOl8  cm-3). 

Fig.  9  Spatial  distribution  of  ion  densities  in  the  incoming  and  reflected  plasma  flows 
at  different  time  moments  (a)  and  ionization  degree  (b)  (one-dimensional  version  n0  =lOl8 
cm-3). 


Fig.  10  Spatial  current  density  distribution  at  some  time  moments  in  the  direction  per¬ 
pendicular  to  the  target. 

Fig.  li  Spatial  distribution  of  longitudinal  component  of  the  generated  magnetic  field. 

Fig.  12  Spatial  distribution  of  ion  temperatures  in  the  incoming  plasma  flow  at  three 
successive  time  moments  (two-dimensional  calculation). 

Fig.  13  Spatial  distribution  of  electron  temperatures  at  three  successive  time  moments 
(two-dimensional  calculation). 

Fig.  14  Spatial  distribution  of  ion  temperatures  in  the  reflected  plasma  flow  at  three 
successive  time  moments  (two-dimensional  calculation). 

Fig.  15  Spatial  distribution  of  ion  densities  in  the  incoming  plasma  flow  at  three  suc¬ 
cessive  time  moments  (two-dimensional  calculation). 

Fig.  16  Spatial  distribution  of  ion  densities  in  the  reflected  plasma  flow  at  three  suc¬ 
cessive  time  moments  (two-dimensional  calculation). 

Fig.  17  Field  of  ion  velocities  in  the  reflected  plasma  flow  at  three  successive  time  mo¬ 
ments  (two-dimensional  calculation). 

Fig.  18  Field  of  electron  velocities  at  three  successive  time  moments  (two-dimensional 
calculation). 


User’s  guide 


Initial  program  text  is  in  file  BUNCH.FOR  (programming  language  is 
«FORTRAN-77»).  The  initial  data  are  introduced  from  file  OTL.PAR,  where  ANX  is  the 
characteristic  meaning  of  density  (cm-3),  Tioo,  T200,  TEoo  are  the  initial  values  of 
ion  temperatures  in  the  incoming  flow,  reflected  flow,  electrons  (eV),  U10,  U20  are  the 
initial  velocity  values  (in  dimensionless  form),  AA  is  the  initial  clot  radius  (cm) 
(RX=AA),  VX  is  the  characteristic  velocity  value  (cm/s). 

In  the  program  text  itself:  NX,  NR  are  the  number  of  cells  in  the  calculation  re¬ 
gion  of  x  and  r,  MSC  determines  the  number  of  time  steps  for  data  output  to  files: 
*.DAT,  MMC  determines  the  number  of  time  steps  for  control  data  output  to 
BUNCH.LOG  or  on  screen  (  1st  column  is  the  number  of  time  steps,  2nd  is  the  dimensionless 
time  (tx=Rx/Vx),  3rd  is  the  dimensionless  time  step). 

The  data  output  to  files  *.DAT  is  in  the  form  (for  example): 

WRITE(N,*)  T  is  the  time  (s), 

WRITE(N,*)  NXX,NRR  is  the  number  of  cells  in  x  and  r  (pa3Mepbi  MaccuBOB 
nepeMeHHbix), 

WRITE(N,*)  (XX(I),I=o,NX)  are  the  coordinates  of  cell  boundaries  on  x  (cm), 
WRITE(N,*)  (RR(J),J=o,NR)  are  the  coordinates  of  cell  boundaries  on  r  (cm), 
WRITE(N,*)  ((Ti(I,J),I=o,NX),  J=o,NR)  is  the  array  ion  temperature  values  in  the 
incoming  plasma  flow. 

Output  characteristics  (data  arrays): 

Ti(I,J)  is  the  ion  temperature  in  the  incoming  flow,  eV 
T2(I,J)  is  the  ion  temperature  in  the  reflected  flow,  eV 
TE(I,J)  is  the  electron  temperature,  eV, 

Ui(I,J)  is  the  Ux  component  of  the  vector  of  ion  velocity  in  the  incoming  flow,  *io8 
cm/c, 

Vi(I,J)  is  the  Ur  component  of  the  vector  of  ion  velocity  in  the  incoming  flow,  *io8 
cm/c, 

U2(I,J)  is  the  Ux  component  of  the  vector  of  ion  velocity  in  the  reflected  flow,  *io8 
cm/c, 

V2(I,J)  is  the  Ur  component  of  the  vector  of  ion  velocity  in  the  reflected  flow,  *108 
cm/c, 


UE(I,J)  is  the  Ux  component  of  the  vector  of  electron  velocity,  *io8  cm/c, 

VE(I,J)  is  the  Ur  component  of  the  vector  of  electron  velocity,  *io8  cm/c, 

ANi(I,J)  is  the  ion  density  in  the  incoming  flow  *(nx=io17  cnr3), 

AN2(I,J)  is  the  ion  density  in  the  reflected  flow  *(nx=io17  cm-3), 

AALFS(I,  J)  is  the  ionization  degree  a, 

LIN(I,J,IZ,M,K)  is  the  array  of  radiation  coefficients  Ezmk,  in  cells  (i,j)  for  ions  with 
charge  Z  in  transitions  m-»k. 

To  plot  the  isoline  graphs  the  program  IZOLD.FOR  is  used.  The  isoline  values 
may  be  taken  automatically  (uniform  decomposition  between  minimum  and  maxi¬ 
mum  meanings  in  data  files),  or  by  handle  (on  users’s  desire).  To  do  this,  the  file 
GRAFCON1-3.DAT  is  used  (three  graphs  are  plotted  on  sheet  (see  the  examples  of 
twoodimensional  calculations)),  7th  line  of  this  file  takes  the  form:  1  (or  o  if  automati¬ 
cally  mode  is  used),  7  number  of  isolines,  50  100  200  300  500  750  1000  -  (  meanings, 
for  example  Tx,  for  isoline  with  a  proper  number). 

To  plot  a  field  of  velocities,  the  program  VEL.FOR  is  used.  During  its  opera¬ 
tion  first  the  file  name  should  be  introduced  in  which  the  meanings  of  1st  coordinate 
Ux  of  velocity  vector  are  placed  (for  example,  Ui.DAT),  and  then  the  file  name  with 
meanings  of  2nd  coordinate  Ur  of  vector  (for  example,  Vi.DAT).  The  file 
GRAFCON.DAT  is  necessary  for  program  operation. 

Following  translation  of  files  IZOLD.FOR  h  VEL.FOR,  it  is  necessary  to  join 
the  program  library  GRAFOR  by  the  command: 

$  LINK  IZOLD,SYS$LIBRARY:GRAFOR/LIB. 

Program  translation,  calculations,  graph  plotting  were  carried  out  on  mini¬ 
computer  with  ALPHA  processor  and  with  operation  system  VAX  VMS. 
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